Time-Reversal Coherent Control in Nanoplasmonics 
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We introduce an approach to determining the required waveforms to coherently control the optical 
energy localization in plasmonic nanosystems. This approach is based on the impulsive localized ex- 
citation of the nanosystem and time reversal of the generated far-zone field at a single point with one 
polarization. Despite strong interaction and significant dephasing and dissipation in metal plasmonic 
systems, and incompleteness of this time reversal, the proposed approach proves to be very efficient 
in controlling the nanoscale optical fields. Possible applications include nanoscale spectroscopy and 
photomodification, ultradense memory, and information processing on the nanoscale. 
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One of the challenging problems is control of the 
nanoscale localization of optical excitation energy in 
nanosystems. The principal physical origin of this prob- 
lem is that the wavelength of the optical radiation is or- 
ders of magnitude larger than the size of a nanoscale sys- 
tem. Therefore the optical radiation cannot be focused 
into a nanoscale spot, i.e., it does not posses spatial de- 
grees of freedom on the nanoscale. However, it does pos- 
sess the temporal degrees of freedom, which is the same, 
in a different language, as the frequency or phase de- 
grees of freedom. The idea to use the phase modulation 
of the optical radiation to control the nanoscale localiza- 
tion of the optical energy has been proposed^ and led to a 
significant subsequent development both theoretical 2 -^'^ 
and experimental^. In symmetric systems, polarization 
of the excitation has been an efficient degree of freedom 
complementing the phase and amplitude modulations^ 7 - 
Coherent control of quantum and nanoscopic systems is 
a powerful tool of defining pathways of optically-excited 
processes in them^ i 9 ' 10 : 11 ] 12 : 13 

One of the fundamental problems in the coherent con- 
trol is the solution of its "inverse problem" : finding an 
optical waveform that sends the controlled system along 
the required excitation pathway. One approach to this 
problem is the adaptive optimum control that has proved 
successful in a wide class of problems Sil i 10 : 11 However, in 
the adaptive algorithms, it is sometimes difficult to inter- 
pret the obtained complicated waveforms. A problem in 
theoretical investigations is that implementations of the 
adaptive algorithms are often computationally costly. 

In this Letter we propose and theoretically investigate 
a novel approach to finding an efficient optical pulse con- 
trolling a nanosystem. It is based on an idea of time 
reversal (or, back propagation). We start with an initial 
state of the plasmonic nanosystem where a localized ex- 
citation is prepared at a desired nanosite. Then we solve 
the direct problem of evolution and propagate the fields 
to the far zone. At some moment t r , we time- reverse the 
far-zone pulse and send it back to the system. If the sys- 
tem were completely time reversible, then its evolution 
would back-track itself causing the concentration of en- 
ergy at the desired site in time t r after the instance of 



the time reversal. This idea is significantly based on the 
previous acoustic and microwave studies i 14 d 5 i 16 

However, there are principal differences of surface 
plasmon (SP) eigenmodes in metal nanosystems from 
the reverberating, leaky modes in acoustics and mi- 
crowaves of Refs. Hflalla . which make controlling the 
SPs much more difficult. First, due to the strong in- 
teraction of dipolar excitations on the nanoscale, the 
SPs in nanosystems form chaotic eigenmodes that can 
be delocalized over the entire nanosystenii 17 : 18 This phe- 
nomenon manifests itself as correlated "hot spots" of lo- 
cal fields. Obviously, both the delocalization of the SPs 
and their chaoticity (including high sensitivity to param- 
eters) hamper the ability to control them. Second im- 
portant difference is that the real metals in the plas- 
monic spectral region are lossyS^ Therefore, SPs have 
finite life times and are not exactly time-reversible. The 
third problem, which is common to the acoustics, mi- 
crowaves, and plasmonics, is that the far field does not 
contain the full information of the internal state of the 
systems: the evanescent fields are vanishingly small and 
lost in the far zone. Moreover, the far-zone field is nor- 
mally measured in a limited number of points with in- 
complete polarization information (in the extreme case, 
at one point with a single polarization). Therefore, the 
full time-reversal of the field is impossible. These serious 
problems notwithstanding, as we show below in this Let- 
ter, the time reversal of the far-zone field even at a single 
point with a single polarization produces a signal that is 
capable of providing an excellent control of the optical 
field nano-localization. 

Turning to the theory, we consider a nanostructured 
system consisting of metal and dielectric with the per- 
mittivities e m and Ed, respectively. The entire size of 
this system is assumed to be much less than the wave- 
length of the excitation radiation. Therefore we can use 
the quasistatic spectral theoryS^ We start with an op- 
tical dipole d(ro,t) localized at a point ro at the metal 
surface whose density is P(r) = 6(r — ro)d(ro, t), and the 
dependence on time t is a short pulse containing frequen- 
cies ui centered around the carrier frequency ojq. 

This initial oscillating dipole causes the appearance of 
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FIG. 1: (a) Geometry of the random planar composite used 
in the computations, shown in the xz projection. The unit 
length for the axes is 1 nm, but can be scaled. The system is 
scalable in the limits allowed by the quasistatic approxima- 
tion, (b) Average (over pulse time) intensity of local fields 
is displayed distributed over the surface of the nanosystem. 
The fields magnitude is shown in the units of the excitation 
pulse amplitude, which is set as 1, and whose length is 5 fs 
with mean frequency ujq = 1-2 eV. 



local fields E L (r, t) in the system that are given by^i 

t 47T 

££(r,o/) = — G^(r,r ;o; )dp(r ,w) , (1) 
£d 

where the Greek subscripts denote vector indices with 
summation over repeated indices implied. Here and be- 
low, by indicating a frequency variable lo we imply the 
Fourier transform of the corresponding temporal func- 
tion; e.g., E(r, u) = E(r, t) exp(iu>t)dt. Retarded 
dyadic Green's function G r a p is expressed in terms of the 
corresponding scalar Green's function G r : 



dr a dr' 



r G r (r,r';c) . 



(2) 







This is given as an expansion over the eigenfunctions ip n 
and eigenvalues s n of the SP eigenproblem ^ 20 ! 21 



<r(r,r»=£ 



<p n (r) (f n (r')* 
s(lo) - S n 



(3) 

The total radiating dipole moment of the nanosystem 
D, which defines the field in the far zone, is the seed 
dipole d plus the dipole of the entire system induced 
by field E L ((T|) that renormalizes and enhances it (the 
antenna effect). It is given in the frequency domain as 21 



D a (to) 



Sal: 



1 



s{uj) 



(4) 



(? a/3 (r,a,) = / ^(r,r'; U )0(r')dV , (5) 
Jv 

where 6(r') is the characteristic function equal to 1 when 
r' belongs to the metal and to otherwise. To complete 
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FIG. 2: (a) Excitation field at the system as computed by 
time-reversal for the initial dipole at the A point (normalized 
to 1 at the maximum), (b)-(d) Distributions of the local field 
intensity at the points A, B, and C, correspondingly, calcu- 
lated for the instances of their respective maxima, (e) Time 
evolution of the local intensity at the targeted point A. (f) 
Distribution of the time-averaged squared intensity over the 
surface of the nanosystem. Units of field and intensity are 
arbitrary but consistent for all the panels. 



the solution of the inverse problem of the coherent con- 
trol, we find the time-reversed radiating dipole moment 



D*{uj) exp(— iu>t)dui/ (2w) 



(6) 



The field generated by this oscillating dipole in the far 
zone is then used to excite the nanosystem. 

The final step is the solution of the direct problem, i.e. 
finding the the local field E(r, t) in the system that is 
excited by the time-reversed uniform field ~E T (t). This 
can be expressed in the Fourier domain as 



E a (r, uj) = [S al 3 + g a p(r, w)] E^{uj) 



(7) 



If the time-reversal is efficient in solving the inverse prob- 
lem, then the local field, Fourier-transformed to the time 
domain, should demonstrate the concentration at the ini- 
tial site ro in time t r after the reversal. 

As a numerical illustration, we consider a random pla- 
nar composite which is a 4-nm thick layer of silver— in 
vacuum. This layer is 50 percent randomly filled with 
2x2x2 nm 3 unit cells, as shown in Fig. Q] (a). The 



total size of this composites is 32 x 4 x 32 nm 3 . In the 
quasistatic approximation, the system is scalable, as long 
as its total size is still much less the light wavelength. 

There are a priori limitations on the coherent control 
in the plasmonic nanosystems. In particular, the local 
field energy can only be localized at the sites where the 
eigenmodes with frequencies within the bandwidth of the 
excitation pulse are localized. To get an idea where in 
the nanosystem such a localization takes place, we apply 
a very short, 5-fs duration, unmodulated Gaussian exci- 
tation pulse whose carrier frequency ujq = 1.2 eV is in the 
window of the least dephasing of the SPs^ The result- 
ing local field intensity averaged over the pulse time, is 
shown in Fig.[T](b). Among the many hot spots of the lo- 
cal fields seen in this panel, we choose three peaks marked 
as A, B, and C, whose p = [x, z) coordinates at the sur- 
face are p A = (8,20), p B = (5,5), and p c = (20,28), 
correspondingly. In all cases the time dependence of the 
initial dipole d(ro,t) has been set as a pulse with a 20 
fs Gaussian envelope and ujq = 1.2 eV. The initial dipole 
was polarized in the z-direction. Separate computations 
for this dipole x-polarized have given very similar results 
(not shown), which is due to the strong depolarization ef- 
fect of the present random nanostructure. This property 
of random, complex nanosystems is in sharp contrast to 
the polarization-driven control for a symmetric system.— 
The time-reversed excitation field in the far zone has been 
calculated from the dipole ([6]). Its maximum amplitude 
on the system has been normalized to 1 to make the com- 
parison of the plasmonic enhancements easier. 

When the initial dipole is at the point A, the calcu- 
lated time-reversed electric field incident on the system 
is shown in Fig. [2] (a) where t r — 536 fs. We can see 
that this field is dramatically different from the seed 20- 
fs Gaussian-envelope polarization that generated it (af- 
ter the time reversal). The pulse is relatively long, with 
bursts of fields and their revivals, accompanied by the 
general decay due to the dephasing. This pulse looks 
very similar to the pulses obtained by the time reversal 
in acoustics and microwave o 14 i 15 ' 16 where this behavior is 
due the reverberations of the propagating waves repeat- 
edly reflected from the boundaries and inhomogeneities 
of the system; the decay of the signal is due to the leak- 
age of the wave energy from those open systems. How- 
ever, the similarity stops here, because in our case the 
observed beatings are due to the interference of the lo- 
calized, non-propagating (quasistatic) SPs eigenmodes; 
the decay is due to their dephasing occurring both due 
to the multitude and randomness of the eigenmode fre- 
quencies involved, and also due to the dephasing of the 
metal electron polarization as described by Ime m . 

In Figs. (b)-(d) we display the local field intensity 
I{r,t) = |E(r,i)| 2 distributed over the surface of the 
nanosystem for three moments of time where the intensi- 
ties at sites A, B, and C are maximum, correspondingly. 
Note that in all these cases the intial dipole is at the A 
point. The maximum concentration of energy at the A 
site is almost perfect [panel (b)]; it is reached at t = 536 
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FIG. 3: The same as in Fig. [2] but for the initial dipole at the 
B point. 



fs, which coincides with the expected time t r (the end of 
the excitation pulse). Comparing to the case of an um- 
modulated pulse [Fig. [1] (b)], the excitation of the other 
peaks is almost completely suppressed. Even when the 
"undesired" peaks B and C go through their temporal 
maxima [panels (c) and (d)], the targeted peak A is still 
dominant. Not only the spatial structure of the local 
fields is highly concentrated. Also the temporal evolution 
of the local field intensity at the A point shown in Fig. 
(e) is restored almost completely to its initial Gaussian 
envelope (though with some pedestal). Finally, we dis- 
play in Fig. [2] (f ) the time-integrated square of the inten- 
sity, which describes the distribution of the two-photon 
electron-emission current as measured, e. g., by a photoc- 
mission electron microscope (PEEM). 6 ' 7 Such a current 
is almost ideally concentrated at the targeted point A. 

Similar case for the initial dipole at the B point is il- 
lustrated in Fig. [3) The excitation, time-reversed pulse 
[panel (a)] is significantly different from the previous case: 
the beatings are obviously much less frequent. The spa- 
tial distributions of the local field intensity displayed in 
panels (b)-(d) reveal that the initially weak peak B [cf. 
Fig. [T] (b)] is relatively very much enhanced. It reaches 
its maximum [panel (c)] at t = 536 fs coinciding with the 
expected (back-tracking) time £ r , where it is the largest 
peak]. This is certainly a success of the coherent control. 
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FIG. 4: The same as in Fig. [2] but for the initial dipole at the 
C point. 



Moreover, the time evolution of the local fields at the B 
site [panel (e)] shows an excellent temporal concentration 
and reproduction of the seed 20 fs Gaussian pulse. How- 
ever, as panel (d) shows, the maximum magnitude of the 
peak at the C site is comparable, though smaller, than 
that of the maximum B peak. It is also important that 
the initially strongest peak A [Fig. Q] (b)] is significantly 
suppressed. The time-averaged nonlinear current [panel 
(f)] is nevertheless dominated by the C site, which is due 
to the very long-lived local fields at that site. Thus the 
temporal concentration at the targeted point B driven by 
the time-reversed field is sharp and transient in time. 



For the C point as targeted, illustrated in Fig. [4] the 
excitation, time-reversed pulse [panel (a)] is very long, 
lacking strong beatings, which implies a weak dephasing. 
The C peak certainly dominates the temporal dynamics 
reaching its maximum panel (d)] at t = 533 fs, which 
is just one period of oscillations shifted from the back- 
tracking time t r — 556 fs. The C peak also dominates the 
nonlinear current [panel (f)]. The temporal dynamics at 
the targeted C point [panel (e)] shows a rather broadened 
peak, but its very center exhibits some narrow spike at 
t « t r . Overall, the time-reversal coherent control is very 
efficient at the selective concentration of the excitation 
energy at this point. 

To briefly conclude, we have introduced an efficient 
approach to solving an and important and formidable 
problem of the coherent control of the local optical en- 
ergy distribution in plasmonic nanosystems. We have 
start with a localized dipole producing a short pulse of 
optical oscillations at a targeted site of the nanosystem. 
The field of this dipole with one polarization at a sin- 
gle point in the far zone is time-reversed and used as an 
excitation pulse. We have shown above that despite the 
significant problems in time reversing a lossy, strongly 
interacting plasmonic nanosystems using incomplete in- 
formation, it is still possible to efficiently concentrate the 
energy of ultrafast optical fields at the targeted nano-site. 
Thus, the time reversal provides a powerful method to 
solve this fundamental problem of the coherent control. 
This method can be used either as an alternative to or 
in combination with the adaptive coherent control where 
it provides the initial pulse. This proposed approach can 
be used for controlling the ultrafast local optical dynam- 
ics in nanosystems for a variety of applications, including 
superdense and ultrafast optical memory and computing 
on the nanoscale, ultrafast local spectroscopy and photo- 
chemistry on the nanoscale, and others. 
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